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The Unipolar Induction 

P. Hrasktfl 

1 The phenomenon of unipolar induction 

Consider a permanent magnet in the form of an electrically conducting disk 
magnetized along its symmetry axis and rotating around it with constant an- 
gular velocity O. If, using sliding contacts, one connects the axis and the rim 
of the disk by means of a rigid linear conductor electric current will go through 
the latter (unipolar generator, Fig. 1). An electromotive force arises even when 
the disk is made of a dielectric not supporting electric current. The mechanism 
which gives rise to this electromotive force is called unipolar inductioiQ 
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Figure 1 



How to calculate the magnitude of the effect? When the disk is at rest and 
the conductor rotates in the opposite direction with angular velocity — J7, it 
crosses contiuously the lines of force of the magnet and, therefore, an electro- 
motive force 

Sad = I [(r x (-«)) x B] ■ dr (1) 

J ABCD 

arises in it through motional induction. Since the relative motion is the same 
in both cases one might suppose that this same electromotive force arises also 
when the disk rotates together with the lines of force attached to it and the 
conductor remains at reslo- 

This reasoning is, however, of highly dubious value for at least two reasons. 

First, rotation is absolute and it is by no means immaterial whether the disk 
or the conductor is rotating. In unipolar induction the rotation of the disk is 
supposed to take place in an inertial frame of reference. When, on the other 
hand, the disk is assumed resting, the point of view is changed and the reference 
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frame becomes rotating. In mechanics a change of this kind is known to be 
accompanied by a corresponding change of the equations of motions, consisting 
in the inclusion of centrifugal and Coriolis forces. There is absolutely no reason 
to expect that Maxwell-equations remain intact under the same transformation 
(See Section 3 below). 

Second, the lines of force are of a purely mathematical device, they motion 
can be neither observed nor calculated. It is, therefore, meaningless to differen- 
tiate, by reference to the motion of the lines of force, between the magnetic fields 
of magnets which rotate with different angular velocity around their symmetry 
axis. Hence in the two cases compared above the magnetic field is actually the 
same and the linear conductor which was originally at rest is simply replaced by 
a rotating one. These two situations are, however, completely unrelated since 
they differ from each other more than in a mere change of the point of view. 

The correct explanation of the phenomenon is based on the observation, 
that the elements of the magnet rotate in the field of all the other elements. 
The magnet, therefore, rotates in its own magnetic field and — due to motional 
induction — an electromotive force 

(V x B) = [(r x n) x B] 

will act in the element of it which at the given moment is in the position r. Since 
the magnet is assumed conducting, the free charges get moved by this force until 
the field E due to the volume and surface polarization charges compensates the 
induced electromotive force within the magnet: 

E = -(VxB). (2) 

It is the electric field of these polarization charges outside the magnet which 
drives the current in the contour ABCD. Therefore, instead of ([1]), the correct 
electromotive force is given by the formula 

Sab = f E dr. (3) 

J ABCD 

The electric field, surrounding the magnet, remains obscured in the previous 
(erroneous) explanation^. But, in spite of their fundamental difference, both 
the correct and the erroneous approaches lead to exactly the same electromotive 
force. This follows from the general forrrj^l 

* = - I [E + (V x B)] ■ dr (4) 

Jc't 

of the law of induction which is valid even for moving contours C*. 

4 Since the polarization charges rotate together with the magnet they generate a constant 
magnetic field which modifies the original field of the magnet. It is easily seen that this 
secondary field is proportional to Q 2 . In what follows we will confine ourselves to first order 
in Q in which this modification of the magnetic field plays no role. 

5 The proof of this formula is found in Appendix 1. 
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Assume that the contour ABCD rotates together with the magnet and close 
it by adding a part DP A within the magnet. Due to axial symmetry the mag- 
netic flux through this closed contour is constant in time, therefore 



* [E + (V x B)] ■ dr = 0. 

J ABCDPA 



However, as a consequence of ©, the integrand on the part DP A vanishes and 
we have 



L 



E-dr = (V x B)-dr 

ABCD J ABCD 



= - / [(r x n) x B] ■ dr 

J ABCD 

which proves the numerical equality of the two forms of Sab- 

The first method, however, is by no means justified by this coincidence. 
When, for example, the magnet is a dielectric the first method leads obviously 
to the same electromotive force as before while the correct electromotive force 
given by the second method turns out to depend on the dielectric constant of 
the disk (see next section). 

Summing up: Unipolar induction consists in the electric polarization of the 
permanent magnet, rotating around its symmetry axis. In the special case 
of magnets made of conducting material the formula for the potential drop 
along contours like ABCD coincides with the formula for the electromotive 
force induced in a contour rotating around a magnet which is at rest in an 
inertial frame of reference. In a practical calculation however, instead of ([5]), 
one can start from the known inside electric field © and calculate the outside 
field from the continuity of the tangential component of E across the surface 
of the magnet and the absence of the net charge on it. For a spherical magnet 
of radius a and constant magnetic dipole density M this procedure leads to the 
potential 

* = -§aW^M^) (r > a) (6) 
9 v 

outside the sphere. 



2 The role of the motionally induced electric dipole density 



The magnetic dipole density M , moving with velocity V , acquires electric 
dipole density Pm- To linear order in O 

Pm = \(V x M). 

c z 

Usually this phenomenon is treated in the framework of relativity theory though, 
being proportional to V/c 2 rather than (V/c) 2 , it is not a genuine relativistic 
effect. Does this electric polarization contribute to unipolar induction? 

6 See paper The rotating magnet on my homepage www.hrasko.com/peter . 
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It does, but if the magnet conducts electricity, the consideration of the previ- 
ous section remains unaltered. In stationary rotation the electric current density 
J vanishes in both conductors and dielectrics. According to Ohm's law we have 
J = a\E + (V x B)), In conductors a ^ 0, the condition J — leads to 
the electric field ([2]) inside the magnet and thereby the field outside the magnet 
becomes unambiguously determined. It is of no significance whether this field is 
generated solely by the free charges or the polarization charge density — div Pm 
also contributes to it. 

When, on the other hand, the magnet is made of dielectric, in the Ohm's law 
we have a = and J vanishes without requiring ((2]) to be fulfilled. In this case 
the relevant electrostatic problem can be solved 6 , starting from the equation 
div D — where now 

D = eE - e x(V x B) - P M - 

The second term on r.h.s. is the polarization generated by the motionally in- 
duced electromotive force (x is the electric susceptibility). 

We arrive again at a well defined electrostatic problem which is easily solved 
for a spherical magnet. For the potential difference between the axis of rotation 
and the "equator" the formula 

A$ ms = A$ corld • 2g + £ ° 
2e + 3e 

is obtained ($ cond i s given in ©). As e — > oo we have A$'" s — > A$ cond 

as expected. When the magnet is neither conducting nor polarizable we ob- 
3 

tain A$"lf eo = -A$ co " d which is the contribution of P M alone to unipolar 
induction. 

3 Transition to the corotating frame 

Reference frame must always be chosen so as to simplify the analysis of 
the problem under study. Since there are important problems which are most 
conveniently studied in the rest frame of the rotating magnet we turn to the 
description of the unipolar induction in this latter frame. Observable properties 
such as whether an electric bulb in the contour ABCD of Fig. I is or is not 
gleaming cannot, of course, be altered by mere change of the point of view but 
to verify this we must first transcript Maxwell-equations into rotating frame. 

In order to anticipate the necessary modifications assume that in the immedi- 
ate vicinity of the axis of rotation a constant linear charge density of magnitude 
A is concentrated. Then, if the equation div EqE — p remained valid in the 
rotating frame too, the flux of E through the multitude of cylindrical surfaces, 
surrounding the axis of rotation, would remain equal to X/sq per unit height in 
both rotating and inertial frames and this requirement would strongly limit the 
possible change of the electric field in the transition from an inertial frame to 
the corotating one. 
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But such a change must take place. The Coulomb force eE by which the 
charge density A acts on a point charge e is certainly different when the charge 
is at rest with respect to the one or the other frame. This is suggested by the 
transformation law 

E' ^E + (V xB) (V < c) 

between two inertial systems. Though inapplicable when transition takes place 
to rotating frame, this equation suggests that the electric field does change in 
the latter transition also, leading an observable difference in the value of the 
Coulomb force in the two frames. 

Let us assume now that an electric current I flows along the axis of rotation. 
Than, if the equation rot B = (IqJ remained valid in the rotating frame, the field 
B would presumably also have to remain essentially the same. But this would 
almost certainly be in conflict with the relativistic transformation properties of 
the fields and the expected modification of the Lorentz-force. 

The way out of this dilemma is suggested by the electrostatics of polarized 
media: On the level of field equations the fields E, B which determine the 
force exerted on charged particles must be distinguished from the fields D, H 
generated by the latters. At the same time, an appropriate pair of "material 
equations" must be postulated in order to restore the correct number of degrees 
of freedom. The whole process must naturally be based on the known relativistic 
transformation properties of the field in the form of an antisymmetric field tensor 

The shortest way to this form of electrodynamics in a rotating frame goes 
through general relativity where the problem has been solved in complete gen- 
erality for an arbitrary system of coordinate^. Specification of the general 
formulae leads to the following equations valid in rotating coordinate system: 

9B „ n 

divs-o ~dr +rotE = 



divD^p dD_ 

dt 

The "material equations" are 



(7) 

rot H - J. 



D= , £ ° =E+\(HxV) (8) 



B= ; Mo =H +L(VxE). (9) 



V(r) is a vector field present in the rotating frame given by the formula V 
(r x ft). In the rotating frame Lorentz-force acquires an additional term: 



Fl = eE + e(v x B) — %{E ■ v) V . (10) 



7 See L.D. Landau and E.M. Lifshitz The Classical Theory of Fields Fourth Revised English 
Edition, § 90, p. 275, and Appendix 2 at the end of the present work. 
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Apply now these equations to unipolar induction. To this end we have to 
replace H in ([9]) by the sum H + M within the magnet. The jump conditions 
on the surface of the magnet follow from ([7]) and are the usual ones: Normal 
components of D and B and tangential components of E and H are continuous 
which implies the continuity of the electric potential. 

The same argument which in Section 1 led to the relation @ now gives 
E = inside the magnet. Outside we have divD = p = which leads to the 
equation 

= -\div(H x V). (11) 

This formula is the starting point of the calculation of the electric field in the 
corotating frame. 

As we have mentioned earlier 4 our considerations in the inertial frame were 
implicitely restricted to linear order in f2 since back reaction of the induced 
currents on the magnetic field was neglected. In the rotating frame no such 
currents arise and nonlinearity in Q is made explicite in the equations f7|-(|9|) 
which is a definite advantage of the rotating frame. 

In order to confine ourselves to linear order agaird we note that when Q = 
the electric field E obviously disappears and so it is proportional to f2. Since 
the same is true for the field V(r) the second term on the r.h.s. of |9]) is of the 
order of f2 2 and must be neglected. At the same time the square roots in the 
"material equations" and in (fTTj) must be replaced by 1. 

In this approximation, therefore, the magnetic field in the corotating system 
is exactly the same well known dipole field as in the laboratory (inertial) fram^|. 
It can be substituted into the r.h.s. of (fTTj) which can than be used to calculate 
E. Since the fields are static we have from (0 rot E = and an electric potential 
<]>' may be introduced into (fTTj) by the usual relation E = — grad<I>': 

A$' = \div(H x V). (12) 

Cr 

This equation must be solved under the following conditions: (1) Since the 
tangential component of the electric field on the surface of the magnet is con- 
tinuous and, therefore, equal to zero $' must be constant along this surface; (2) 
Since the magnet is uncharged the flux of the field D across the closed surfaces, 
surrounding the magnet, must be zero. For the spherical magnet the solution is 

1 sin 2 ?9 

$' = $ _ -a 3 n MCl =— - (r > a) (13) 
3 v 

which is to be compared with ((6]). 



"This limit can be directly contrived without making the detour through general relativity 
(see Landau and Lifshitz Electrodynamics of Contiuous Media § 76). 

9 In the literature on the magnetic hydrodinamics of rotating fluids this fact is often ex- 
pressed as e.g. "rigid body rotation... rotates a magnetic field without distortion" (H. K. 
Moffat Magnetic Field Generation in Electrically Conducting Fluids Cambridge f978, p. 53) 
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Let us calculate finally the electromotive force £' AB in the contour ABCD 
of Section 1 which rotates with angular velocity lu = — Q in the corotating 
system. This force is determined by the Lorentz-force equation (|10[) in which 
now v — (r x u>) = — (r x f2). The third term of the force is of the order of 
r2 2 and must be neglected while the first term's contribution vanishes since the 
endpoints of the contour lie on the surface of the magnet which is equipotential. 
We have, therefore 



which is exactly the same expression as (pj and J5]). When the magnet is 
insulating its surface ceases to be equipotential and the first term of (fTU| also 
contributes to the electromotive force. This contribution is determined directly 
by fi, being independent of v. 

4 Application of the theory to the Earth 

Erth's core is a huge rotating magnet and it is tempting to apply the ideas 
of the preceding sections to it. To start with let us assume that (1) the core is a 
permanent magnet with nonzero conductivitj0, (2) the magnetic field outside 
the core is a pure dipole field, (3) the geographical and magnetic poles coincide, 
(4) the mantle is electrically neutral (/i = /io, e = Eq, <j = 00, and finally (5) 
the atmosphere is absent. 

These are mostly very crude assumptions but they permit us to apply the 
formulas of the previous sections without further refinements. The electric field 
at a given point on the surface of the Earth can be calculated from (11311 . As- 
sumption (2) allows to replace the unknown polarization density M with the 
empirical value B of the magnetic field on the Equator, equal to B — 33 [iT. On 
the Equator the electric field given by (|13p is a purely radial one and is given 
by the formula 



where R and R c are the Earth's and the core's radiuses respectively. Substitut- 
ing R c w R/2n 3.10 6 m, = 2tt/24 ft" 1 = 7.3 -lO^s" 1 and B = 33.1CT 6 T wc 
obtain E w 10 mV/m which is in itself of a measurable magnitude but stronger 
fields of atmospheric origin make its observation rather problematic. 

In march of 1996 a spherical 1.6-meter diameter satellite was released out 
from the payload bay of Space Shuttle Columbia during its orbiting at a height 
about h=90 km above the Earth. Its tether a long conducting cable (up to the 

10 The Glatzmaier-Roberts model of the geomagnetic field (G.A.Glatzmaier and P.H.Roberts 
Contemporary Physics 38, 269 (1997)) would certainly make possible a much more realistic 
treatment. 

^Electric polarizability of the mantle makes the electric field larger. When the radius of the 
core is half of that of the Earth, the correction factor is equal to (3301e, — 2321)/(670e r +310) 
where e r is the global relative electric permittivity of the mantle. 
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length of 19.1 km) of resistance R, served (among others) to measure the radial 
electric field (called ambient field) in the vicinity of the orbit. The system is 
conveniently called Tethered Satellite System (TSS). 

Free electrons in the ionosphere where TSS operated were attracted to the 
satellite. The electrons travelled along the tether to the orbiter, producing the 
current /. The electric circuit was closed by means of an electron generator 
on the orbiter which returned charged particles back into the ionosphere. The 
product IR gave the electromotive force along the tether the dominant part of 
which was due to the motional induction in the magnetic field of the Eartrr^l 
and could be calculated with sufficiently high accuracy. From the rest of the 
electromotive force the ambient field could be estimated. 
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Figure 2 



Results of the measurement!^! are summarized on Fig. 2. The solid line is 
a theoretical curve calculated from the theory of " tidally driven neutral winds" 
which is outside my competence. The agreement with observations is considered 
satisfactory. 

The dashed curve added by myself has been calculated according to the 
formula 

R 

E r (t) = -2(QRc) P 2 (-sina • sinwi) • B 
R 

derived from Q ; a denotes the angle between the orbital and equatorial planes 
and oj is the angular velocity of the satellite system (the radius R + h of the 
orbit was replaced by R). In the actual experiment their values were a = 25° 
and uj = 27r/90 min^ 1 . The azimuthal angle tot in the orbital plane is measured 
from the intersection with the equatorial plane. 

Though the order of magnitude agreement is quite impressive it must not be 
taken at face value. The curve was calculated under the assumption (5) of no 
atmosphere while the measurement itself was made possible exclusively thanks 
to the existence of the ionosphere. What the order of magnitude of the electric 

12 Since this electromotive force is of outward direction (opposite to the Earth) the orbiter 
had to revolve nearer to the Earth. 

13 S. D. Williams et al. Geophys. Res. Lett., 25 445 (1998). The lower curve shows TSS 
latitude with thick curve indicating local night. 
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field does probably indicate is that this field should be taken into consideration 
among the input data when complicated processes in the ionosphere are studied. 

One may wonder whether unipolar induction in the Earth's core might not 
have some unexpected influence on satellites not specifically designed to observe 
the ambient electric field. Flyby anomalj0 is a phenomenon which, in spite of 
the efforts of its discoverers, remained so far unexplained. In this situation even 
the most improbable explanations deserve attention. A suggestion of this kind 
might be an effect due to the electric field of the Earth. If satellites can somehow 
(perhaps by photoeffect) get electrically charged then the quadrupole character 
of this field might account for the latitude dependence of the observed anomaly. 

Appendix 1: Derivation of (|4])E1 



Let us start with formula 



#(i) = / B • <2£ = <f 
Js Jc, 



B • effi = d> A-dr (14) 
'£ Jc t 

for the magnetic flux \& through the closed material contour Ct , following from 
the definition B — rot A of the vector potential. The contour Ct may vary its 
shape and length in time and is given in parametric form as 

r = r(X,t) (0<A<1), 
r(l,t) = r(0,t), 

where it is understood that, during the deformations, the same value of the 
parameter A remains attached to the same material pointT^l. 

Since the contour is closed the vector function r(X,t) is a periodic function 
of A of period 1: r(A + 1, t) = r(A, t). The velocity of the point at A is equal to 
dr(X,t) 



v = r — 



dt 



Using this representation of the contour, the integral in JT 
in the form 

dxi{\,t) 



can be written 



*(*) = r 

Jo 



A, 



dX, 



where the notation x±, X2, £3 = x, y, z is introduced and summation over re- 
peated indices is understood. In the integrand the components Ai = Ai(r(X,t), t) 
are functions of both t and A. 

For the time derivative of the flux we obtain 



d^ 
~dt 



^= ^ A 







dAj 
dt 



d 

u dt 

dxj 
dX ' 



A 



dxi 
'dX 
d_ 



dX = 

dxj 
dX 



(15) 



dX. 



14 J.D.Anderson et al, Phys. Rev. Lett. 100, 091102 (2008). 
15 H. K. Moffat op. cit. p. 32 

le Since the integral on the r.h.s. of l|4} is not invariant under time dependent reparametriza- 
tion of the contour it is denned only if the paramctrization is given in physical terms (up to 
an arbitrary choice at some given moment of time) . 
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Here 



dAj 
dt 



dAj 
dt 



dAi dx-j 



dA 
dt 



' dxj 



dxj dt wu xjjuj 

is the "material derivative" of the vector potential on the contour. 
In the second term of the integrand 





fdxA 




( dx{\ 


dvi 


dt 


KdX 


= — 
dX 


V dt ) 


~ dX 



hence 

_ r 1 fdA % dx, dAi d Xi 
~dt~J \ dT~dX +Vj 'dx~~dX' 
The third term can be written as 

a dAi 

^(A iVi ) - — Vi , 
the first of which does not contribute: 



A 



dvj 
dX 



dX. 



— (AiVi) dX = (A 1 w 4 )a=i 



and so in the integrand of (|T7)) we have 



{AiVi)x=o = 0, 



dvi 



dAi 



dAj dxi 



1 dx ~ 1 dx ~ 1 d Xl dx 

The formula (|17|1 can then be transformd into 



~dt 



dA 
dt 



dAj dAj 
dxj dx 1 



dxi 



In the integrand 



therefore 



dA, dAj 



dxj 



~dt 



dx, 

dA 
~dt 



— (v x rot A) 



• dr. 



Here 



rot A = B 



and -— — ~E — grad< 
dt 



(16) 



(17) 



(18) 



(19) 



where <j> is the scalar potential. Substituting these into (fT9|) the scalar potential 
gives no contribution and we arrive at the formula ([?]). 

This formula, therefore, is an exact consequence of the Maxwell-equations 
alone but it becomes especially useful when, taking into consideration the Lorentz- 
force too, we interprete the integral in it as the total electromotive force, acting 
in the contour. This interpretation, however, is valid only in quasistationary ap- 
proximation i. e. when the perturbances may be assumed to cross the domain 
of the contour instantaneously. 
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Appendix 2: Maxwell-equations in general coordinates 

Here we summarize the formulae of Landau and Lifshitz 7 , using SI. 

The space-time coordinate system is characterized by the metric tensor gij 
(latin indices take on values 0, 1, 2, 3). The three dimensional metric is denoted 
by 7 Q/ 3 (greek indices are of value 1, 2, 3) which in terms of g^ is given by the 
formula 

lap = ~g a fi + hg a gp, (20) 

where 

i 90a 

n — Soo, g a — • 

5oo 

Further quantities which will be needed are 

7 = det 7 Q/3 g = dct g {j = -hj 

7 a/3 = ~s Q/3 9 a = i a0 g P = -g Oa - 

We have further 

gi k g kj =s{, la(Tl ° = € 

The Maxwell-equations are 

div£? = rot£ = --LA(^B) (21) 

di Vj D = rot if = -LA^yDJ + J, (22) 



in which rot and div are understood in form valid in general coordinates: 

(23) 

The "material equations" are 



(rota}" - — e a & (^L- ^] 



D = ^=E+-(Hxg) 

Vh c 

Vh c 



(24) 



The continuity equation: 



-^^(V7P) + divJ = 0, (25) 



and, finally, the Lorentz-force: 



F L = eE + e(v x B) (E ■ v) g. 

c 
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The connection between 3D and 4D quantities: 



E a — cFqo 



1 



2^7 



01 
a/3 



D a = -e a cVhF 0a f i H a = -^LJLe aPl F 

From these the components E a , B a , D a and H a are obtained with the aid of 
the 3-dimensional metric. 

Specialization to rotating coordinates: In cylindrical coordinates the tranfor- 
mation formulae from the (primed) inertial system to the (unprimed) rotating 
system are 

t' =t, r' = r, i?' = ip' = <p+ Qt 



which leads to the line element 



ds 2 



1 



£!r sin $ 



o , 9 Qr 2 sin 2 i? „ . 
c 2 dt 2 - 2 dip d(ct)- 



-dr 2 - r 2 {dd 2 + sin 2 tfdip 2 ). 
From this line element the metric tensor (and h) can be read off and we have 

fir 2 sin 2 $ 



9t = g# = g v 



It 



1. 



r 2 sin 2 # 



Since the gamma-tensor is diagonal we have 7 QQ = 1/7 QQ and, therefore 

O 



<f = 9* = 0, 



These formulae lead immediately to the equations ([7]) - (flQ|) in which the 
space coordinates can be chosen arbitrarily. 



